Spin relaxation in semiconductor quantum dots 
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We have studied the physical processes responsible for the spin -flip in GaAs quantum dots. 
' We have calculated the rates for different mechanisms which are related to spin-orbit coupling and 

\ cause a spin-flip during the inelastic relaxation of an electron in the dot both with and without a 

magnetic field. We have shown that the zero-dimensional character of the problem when electron 
wave functions are localized in all directions leads to freezing out of the most effective spin-flip 
mechanisms related to the absence of the inversion centers in the elementary crystal cell and at the 
heterointerface and, as a result, to unusually low spin-flip rates. 

m 

(S| ; PACS numbers: 73.23.Hk; 73.40. Gk 

Quantum dots (QD) are small conductive regions in a semiconductor that contain a tunable number of electrons. 
The shape and size of quantum dots can be varied by changing the gate voltage. Besides, the electronic states can be 
significantly modified by a magnetic field applied perpendicular to the plane of the dot. It gives a unique opportunity 
to study the properties of the electron quantum states in detail and manipulate with the electrons in these artificial 
atoms in a controllable way (see reviewsEm) . 

As it is knownB, quantum dots are considered as possible candidates for building a quantum computer. The crucial 
point of the idea is the necessity to couple dots coherently and keep coherence on sufficiently long time scales. In this 
respect, there is a great demand in a theoretical estimation of the typical dephasing (especially spin dephasing) time 
of the electron in the QD. This situation can occur, for example, when the ground state of the two electron system 
corresponds to a singlet state 5 = and an excited state when one of the electrons occupies the higher orbital level 
corresponds to the total spin of the system S — 1 . Then an inelastic relaxation to the ground state of the dot has to 
c~j . be accompanied by a spin flip. As a result, the spin degree of freedom must be connected through some mechanisms 
to the external environment. Moreover, the electron-electron interaction which can be quite important in determining 
the energy levels of the system, is not important at all for the spin-flip process itself. Therefore, we can consider the 
spin-flip processes within the one-electron approach. 

We have studied the physical processes responsible for the spin- flip in GaAs quantum dots. The unit cell in this 
^ ■ materiaLhas no inversion center which gives rise to an effective spin-orbit coupling in the electron spectrum. It is 
knowroEl, that such spin-orbit coupling provides the main source of spin-flip both in 3D and 2D cases. Besides, in 
the polar-type crystal there is a strong coupling of electrons to the bosonic environment through the piezo-electric 
interaction with acoustic phonons. Such coupling can be important for the inelastic electron relaxation in the GaAs 
crystal both with and without a spin-flip. 

We have calculated the rates for different spin-orbit related mechanisms which cause a spin-flip during the inelastic 
^\ , relaxation of the electron in the dot both with and without a magnetic field. The rates are obtained as a function 
of energy transfer. It is generally observed that the corresponding spin- flip rates are by several orders of magnitude 
lower than for any (elastic and inelastic) spin-flip processes for free 2D electrons, i.e. in the case when there is no 
confinement in the plane. The reason for that is as follows. First of all, some spin-flip mechanisms that are effective in 
the 3D case are frozen out in 2D case. It is known that in the 3D case the spin-flip mechanism related to the spin-orbit 
admixture of the p-band functions to the conduction band functions (the so-called Yafet-Elliot mechanismta) is possible 
even in materials with a bulk inversion center. However, in the 2D case the corresponding scattering amplitude may 
, involve only the following term responsible for the spin-flip: <r [p • p'] oc a z , where p,p' are the two-dimensional 
L" ' electron momenta prior to and after scattering (the z-axis coincides with the normal to the 2D plane). Since a z 
. £h ! contains only diagonal components, this mechanism is suppressed in the 2D case. Therefore, in the GaAs-like crystal 
the most effective mechanism of the spin-flip iftrjhe 2D case is related to broken inversion symmetry either in the 
elementary crystal cell or at the heterointerfaceoQ. This brings about terms in the HamiltoniantJ which are linear in 
the two-dimensional electron momentum and proportional to the first power of the small parameter A/E g <C 1, where 
A is the spin-orbit splitting of the valence band of the bulk GaAs crystal arid E g is the energy gap. This leads to the 
quadratic dependence of the spin-flip rate on this parameter in the 2D caseB. 

In this paper we show that in the case of the quantum dot the zero-dimensional character of the problem leads to 
the additional suppression of the spin-flip rate. Namely, the localized character of the electron wave functions in the 
dot leads to a fairly interesting fact that linear in spin-orbit parameter A terms in the matrix elements are absent in 
true two-dimensional case. As a result, the spin- flip rate is proportional to the power of the parameter A which is 
higher than the second one. The contribution to the spin-flip rate which is quadratic with respect to A exists only if 
we take into account the admixture of the electron wave functions of the higher levels of the size quantization in the 
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z-direction, i.e. the weak deviation from the true 2D motion. This contribution is also small because we consider the 
case of a strong transverse quantization when the lateral electron kinetic energy is much smaller than the distance 
between the quantized levels in the z-direction. 

We start with the following effective- mass Hamiltonian which is derived from the Kane model (see Ref. ||) and 
describes the electron in the conduction band in the presence of a magnetic field B parallel to the z-axis (normal to 
the 2D plane) and an arbitrary confining potential U(r), both vertical and lateral: 

p 2 

H = Ho + Hi + Hph + He-ph]Ha = - — + U (r); 

= T~^2 <3 '[ VC/ ' P] + o /o— ^ tT* 7 ' k + \ Yq& ' £ + \9VBV Z B, (1) 
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Here p = —iti'V + (e/c)A is the 3D electron momentum operator (e > 0), m the effective mass, <r the Pauli 
matrices, m cv the parameter of the Kane model with the dimensionality of mass which determines the "interaction" 
of the conduction and the valence bands with other bands, k x = p x (Py — Vz) ancl other components are obtained by 
cyclic permutation of indices, x,y,z are directed along the main crystallographic axes. This Hamiltonian gives very 
accurate description of the conduction band states in the case E <C E g , where E is the electron energy counted from 
the conduction band edge, and all parameters of this Hamiltonian are very well known for GaAsa. 

The last (Zeeman) term describes the direct interaction of the electron spin with the magnetic field, where g is 
electron g-factor and hb is the Bohr magneton. The first and second terms in Hi describe the spin-orbit interaction. 
The first term is due to the relativistic interaction with the electric field caused by the confinement. An impurity 
potential should be also added to U. The second term in Hi is related to the absence of an inversion centre in the 
elementary cell of the GaAs crystal. In the 3D case this leads to spin splitting of the conduction band proportional 
to the cube of the electron momentum. In the true 2D case averaging of the corresponding term along the motion 
of the electron in the direction perpendicular to the 2D plane gives rise to a spin splitting proportional to the first 
power of the 2D electron momentum!!, this splitting being strongly dependent on the orientation of the vector of the 
normal to the 2D plane with respect to the main crystal axes. Below we shall consider only the case when the normal 
is parallel to the [100] direction which is the most frequently used orientation. 

We have also included in the above Hamiltonian the spin-orbit splitting of the electron spectrum due to the strain 
field produced by acoustic phonons. There <p x = (1/2) {u X y,p y }+ — (l/2){u xz , jv}+, where {, }+ means anticommu- 
tator, ipy, (p z are obtained by cyclic permutations, is the lattice strain tensorQ and Vq is the characteristic velocity 
whose value is well known for GaAsu, V = 8x 10 7 cm/ s. 

We have considered all possible mechanisms of the spin-flip which are related to the spin-orbit coupling. These are: 

1) . The spin-flip event is caused by the spin-orbit relativistic coupling to the electric field (polarization) produced 
by the emitted phonon. This process is described by the term proportional to Vb in Eq. (0) . 

2) . The spin- flip event is caused by the spin-orbit admixture of different spin states. This means that in the 
presence of the spin-orbit term in the Hamiltonian the electron " spin-up" state contains actually a small admixture of 
the "spin-down" state. The electron transition with phonon emission between two states with "opposite" spins leads 
to the spin-flip event. This mechanism is related to the first and second terms in the Hamiltonian Hi and the phonon 
provides only energy conservation here. 

3) . In GaAs the electron g-factor (g = —0.4) differs strongly from the free electron value go = 2 which is due to 
the spin-orbit interaction which mixes valence band and conduction band statesu. This leads to a new mechanism of 
spin-flip in GaAs in the presence of an external magnetic field. Namely, the strain produced by the emitted phonon 
is coupled directly to the electron spin through the variation of the effective g-factor: ~ gHBUiz^iB z . 

It is not clear a-priori which mechanism provides the main source of the spin-flip. We found actually that the 
second mechanism related to the admixture of the different spin states is the most effective one. Below we consider 
these mechanisms one by one. 

1). First mechanism. We investigate only the relaxation of the spin component perpendicular to the 2D plane (S z 
component). Then, using the standard expression for the strain tensor due to the acoustic deformationQ, we get for 
the matrix element describing the electron spin-flip transition between spatial states V^i and \&2 with emission of a 
phonon: 

Mu = tG^) fe e y + % e -](*i I \{(p x + ip y ),exp{iQr)} + | , (2) 

where p is the crystal mass density, e the phonon unit polarization vector, Q = (q, q z ) the phonon wave vector. The 
spin decay rate is given by the Fermi golden rule: 
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r T 4 = ^f / (0(92 + 92) I (*i I ^x + ift/),«P(iQr)}+l^ \ 2 S(hsQ~e), (3) 

where e is the energy transfer. We consider transitions between neighbouring discrete energy levels, so that e = 
TisQ « h /m\ 2 , where A is the dot size in the lateral directions. On the other hand, in practically interesting case of 
transitions between low-lying energy levels, the energy transfer is of the order of electron energy itself. This means 
that perturbation theory (including only one-phonon processes) is valid when e 3> ms 2 , s being the sound velocity. 
From this we obtain the condition q z ^> q w 1/A, i.e. the phonon is emitted almost perpendicular to the 2D plane. 
Then, assuming that q z ZQ <C 1 , where zq is the width of the 2D layer in the z- direction, we can easily calculate the 
integral over q z : 
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where $12 = $iV_<I>2 — ^2^— ^± = ^ 'x + and $1, $2 are the lateral electron wave functions. We give here 
the final expressions in the particular case of elliptic quantum dot with confining frequencies lo x and u> y . For the 
transition n x = 1, n y = =>• n x = 0, n y = (e = hto x ) we have: 

m 3 V 2 e 2 [uJZ 

(5) 

In the case of transition n x — 0, n y = 1 =£- n x = 0, n y = (e = Tiuj y ) in the above formula lo x should be replaced by 
ui y and vice versa . With the use of the Vq value mentioned above, the corresponding relaxation time is found to be 
very long (10 -7- 10 _1 s) for a typical energy transfer of the order of 1 10 K. 

2). Second mechanism. We start with the general expression which gives the phonon-assisted transition rate between 
two states 1 and 2 (which can have either the same or "opposite" spins). We consider here the [100] orientation and 
take into account only the interaction with piezo-phonons which is known to be the most effective in polar crystals 
for a low energy transfer: 

r - - T 2^ / 1 / ^MivM(r)Mr) I 2 5(ns t Q-e), (6) 

where hu is the piezomodulus, ehu = 1.2 x 10 7 eV/ cm for GaAs. As before, the phonon is emitted almost perpen- 
dicular to the 2D electron plane. In the case of the interaction with the piezo-phonon this gives rise to an anisotropy 
effect. Namely, the transition rate strongly depends on the orientation of the normal to the 2D plane with respect 
to the main crystallographic axes. The rate is strongly suppressed when the normal is parallel to the [100] axis and 
this is taken into account in Eq.(^) by the multiplier A t (Q) w 2q 2 /q 2 <C 1 which describes the anisotropy effect for 
a transverse phonon (longitudinal phonons give much smaller contribution). There is no such a small factor for the 
[111] orientation and consequently the transition rates for the [100] and [111] orientations of the heterostructure differ 
by more than an order of magnitude. Proceeding as before, we obtain from Eq.(||): 

T 12 = - 2a *V") a J d 2 r X * l2 I V_ \ 2 Xl2 , (7) 

where X12 = $|$2- For reference we also give the expression for the transition rate without spin-flip which follows 
from Eq.(Q). For the transition n x — 1, n y = =$>■ n x = 0, n y — in an elliptic quantum dot (e = %u> x ) we have: 




Tq= ^(3 + ^), (8) 

Zirpn e 

which for the transfer energy of IK has the value f» 5 x 10 8 s" 1 . Again, in the case of transition n x = 0,n y = 1 => 
n x = 0, n y = (e = huj y ) in the above formula uj x should be replaced by uj v and vice versa . 

In the true 2D case and for the [100] orientation we obtain the following spin-orbit Hamiltonian from Eq.([j]) : 

Hi = ft-* J. + *yPy);P = fa) {2mEg) t mcvEg - ( 9 ) 

Note, that for the other orientation [111] used in practice the spin-orbit Hamiltonian has similar structure (i.e. contains 
the same Pauli matrices), so our conclusions presented below also valid for this case. The constant (3 in Eq.(||) depends 
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on the mean value of the p 2 operator in the state described by the wave function Xo(z) of the lowest quantized level in 
the z- direction and fo* GaAs heterostructures takes values between (1 -=-3) • 10 5 cm/s. In Eq.(||) we have dropped the 
so-called Rashba termta which is believed to be much smaller for GaAs heterostructures (and in any case its presence 
does not influence our conclusions because it has a similar structure). The presence of Tli term leads to a nonzero 
value of the spin- flip transition matrix element. At first sight, the scalar product of the spinors $ corresponding to 
the initial spin-up and the final spin-down states should be proportional to the first power of [3. However, in contrast 
to the extended 2D states, in quantum dots we can actually remove the terms linear in (3 from the Hamiltonian by 
the following spin-dependent transformation: 
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We stress that the boundedness of the electron wave functions is essential for this procedure. Then we obtain 
the Hamiltonian which besides the terms with a unit spin matrix contains only the following spin-dependent term: 
(m(3 2 /h)L z a z , where L z = —ihd/d(p is the angular momentum operator. Therefore, in this approximation the correct 
spin functions are the eigenfunctions of the a z operator and there are no spin-flip processes. The above mentioned 
scalar product of the spinors corresponding to "opposite" spin directions is proportional to the third power of (3. The 
corresponding spin-flip rate ~ T^mfi 2 /Tilo) 3 is very small, T -1 ~ 10 _4 s for the energy transfer of the order of IK. We 
give here the result for the case of circular dot (ui x = ui y = ojq) for the phonon- assisted spin- flip transition between 
first excited and ground states: 
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Note that the term cx a x p x p 2 — & y p y p 2 in the Hamiltonian 7ii (Eq.(|l|)) which is proportional to the third power of 
the in-plane momentum cannot be removed by the above mentioned transformation and gives the contribution to the 
value of X12 . Eq. (|7|) , in the first order of the perturbation theory. This contribution can be comparable to that given 

x = 0, riy = in an elliptic quantum dot (e = Huj x ) is: 



by Eq.(hlj). The result for the transition n x = 1,% = 
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where a — {3/E z , E z = (p 2 )/m. This contribution is of the order of Y §(m$ 2 'hio j E 2 Z ) and becomes comparable to that 
given by Eq.([fl]) at Tioj ~ \J m/3 2 E z which is of the order of several K. 

The structure of the transformed Hamiltonian , i.e. the existence of only the a z matrix, suggests an interesting 
anisotropy effect with respect to the direction of the external magnetic field which can be observed for the spin-flip 
processes. Namely, if the direction of the external magnetic field constitutes some angle 8 with the direction of the 
normal (z-axis) , then the true spin quantization axis will be determined by the magnetic field direction and the above 
mentioned spin-orbit term will lead to spin-flip processes with a rate proportional to /3 4 and angular dependence 
sin 2 6. 

There are, of course, contributions to the spin-flip rate proportional to 1 which are either related to the virtual 
transitions to the higher quantized energy levels in the z-direction or to the presence of an impurity potential which 
leads to the nonseparability of the transverse (z) and longitudinal variables in the Hamiltonian (|l]). However, it can 
be easily checked that both effects give actually the small contributions to the spin-flip rate. For example, in the first 
case the rate ~ To(m(3 2 e 3 /ms 2 E z ) is very small (~ Is -1 ), here E z is the distance between the quantized levels in 
the z-direction. For the second case r^j, ~ To(mp 2 /hcu)(Ui mp /E z ) 2 (zo/r c ) 2 , where hco is the typical distance between 
the levels in a dot, zq is the thickness of the 2D layer in the z-direction, r c is the correlation radius of the donor 
potential and Ui mp is a magnitude of this potential which we take to be not larger than Tilo. This rate is obviously 
smaller than that given by Eq.(12). Finally, in the presence of a magnetic field B directed along the z-axis there is 
also a contribution to the spin- flip rate proportional to f3 2 . It appears if one takes into account the finite Zceman 



splitting in the energy spectrum. We give here the final expression for the case of circular dot (u x 



loq) when 



the solutions of the Tio Hamiltonian are well known. These are Darwin-Fock states which are characterized by two 
quantum numbers: n, I. Then for the transitions n = 0, 1 = ±1, f => n = 0, 1 = 0, j we obtain with the use of Eqs.(R,||): 
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where the energy transfer e(l = ±1) = Hoj + huj c l/2 + g^sB, uj = y/u; 2 + (u^/4) and lu c — eB/mc is the cyclotron 
frequency. Note that this rate is of the order of Y^mf} 2 /T%Lo^)(gpBB /fruJo) 2 and for I = +1 equals 5 ■ 10 +5 s _1 for 
HuiQ ~ IK and B = IT (Zeeman energy w 0.3K). 

Thus, for energy transfers of the order of IK the spin-flip time related to the admixture mechanism cannot be 
shorter than (1(T 5 lCT 6 )s. 

3). Third mechanism. To describe it we need to calculate the following antisymmetric combinations in the presence 
of a deformation 



(14) 



between the states 7,7' = ±1/2 of the conduction band. Energy e c corresponds to the bottom of the conduction 
band, mo is the bare electron mass and p x ,Py are the momentum operators in the absence of the magnetic field. In 
Eq.(pT[) we sum over the valence band states (the heavy and light holes and the split-off band). The Hamiltonian 
H(u) (6x6 matrix) describes the valence bands in the presence of strain: 

H(u)=H(0)+H 1 (u), (15) 

where H (0) is the diagonal (6 x 6) matrix whose elements correspond to the energy position of the top of the valence 
(—Eg) and split-off bands {—(E g + Ai) and Hi(u) is a matrix linear in the strain tensor elements which describes 
the interaction of holes with phononsa. This interaction is described by three deformation potential Constantsa in 
contrast to only one constant for conduction band electrons. Expanding Eq.(H4) with respect to H\(u), we obtain the 
following additional term in the Hamiltonian which should be added to Eq.([l]) and is responsible for the spin- flip in 
the conduction band: 

ov f\ ;f,A f ~\ eB z ™o V d(u xz a x + u yz a y )B z 

H so (u) = *D xy (u)— * — JT — m ^ B _ , (16) 

where 77 = A/ (E g + A) w 0.2 for GaAs and d is one of three deformation constants^ which is of the order of several eV. 
Making use of the Hamiltonian ( |l6| ) and calculating the spin-flip rate along the steps outlined while deriving Eq. (||) , 
we obtain: 
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This equation corresponds again to the spin-flip transitions between the neighbouring orbital levels. The rate given 
by Eq.([l7|) (for hu> c ~ e) has the same order of magnitude as that given by Eq.(§) and is very small. 

The results described above are concerned with inelastic transitions between neighbouring quantized energy levels 
in a dot which corresponds to a relatively large energy transfer. The case of spin-flip processes between the Zeeman 
sublevels of the same orbital level was considered earlier in Ref. |l^ by D.Frenkel. This calculation considers the 
spin-flip processes with phonon emission for electrons in spatially localized states in a quantizing magnetic field (the 
quantum Hall effect regime) and deals with only one spin-flip mechanism related to the spin-orbit interaction with 
the electric field produced by the emitted phonon (the first one in our classification). The analog of his result for 
our case can be easily obtained from Eqs.(||J^), assuming that the energy transfer is much smaller than the typical 
distance Tuv between the levels in a dot: e = g^sB huj which in turn means that the phonon wave length is larger 
than the dot size. We consider here only the case of circular dot with confining frequency ujq and states with n = 
and I = 0, ±1 (the ground state and first two excited states). Then instead of Eq.(||) we have: 

= Vii^sBf (| I I +l)ff , (18) 

This rate is very small (~ 10~ 2 -7- 10 _1 s _1 ) for not too large magnetic field (order of IT). The spin-flip rate related 
to the third mechanism is given by {mw/m) 2 (d/Eg) 2 {(iJLBBf/ps b h A ) and has even smaller value. We found that the 
most effective spin-flip mechanism is that related to the admixture of the different spin states. As it was mentioned 
earlier, in the presence of the Zeeman splitting we have the contribution to a spin-flip rate proportional to (i 2 . Then, 
for example, for the transition between the Zeeman sublevels of the ground state of circular dot with emission of the 
piezo- phonon we obtain: 
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where s; and St are longitudinal and transverse sound velocities and while calculating we have used the anisotropy 
factors for longitudinal and transverse piezo-modes given in Ref. [ll]. The ratio of the rates given by Eq.© and 
Eq.(|l8|) is: (ehi4/mVo/h) 2 (mP 2 ms 2 /(Hujo) 4 ) which for Tiajq ~ IK is of the order of 10 6 . So, the spin-admixture 
mechanism is much more powerful. 

In conclusion, we have calculated the rates for the phonon-assisted spin-flip transitions of an electron in a quantum 
dot for all possible mechanisms which are related to the spin-orbit coupling and provide the connection with the spin 
degree of freedom. We have found unusually low values of the spin-flip rates. We show that it is due to the localized 
character of the electron wave functions in a quantum dot which leads to freezing out of the most effective intrinsic 
spin-flip mechanisms related to the absence of the inversion symmetry in the GaAs-like crystals. 
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